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Abstract. - The effect of quantum fluctuations on the current- volt age characteristics of Joseph- 
son junctions and superconducting nanowires is studied in the underdamped limit. Quantum 
fluctuations induce transitions between a Coulomb-blockade and a supercurrent branch, and can 
significantly modify the shape of current-voltage characteristics in the case of a highly resistive 
environment. Owing to the phase-charge duality, our results can be directly extended to the 
opposite overdamped limit. 



Introduction. — A small-capacitance Josephson 
junction is a quantum system with rich dynamics. The 
two conjugate variables are the superconducting phase dif- 
ference (f) across the junction and the charge Q on its elec- 
trodes. Correspondingly, at low temperatures the behav- 
ior of the junction is determined by the competition be- 
tween the Josephson energy Ej and the charging energy 
Ec = (2e)^/2C, where C is the junction capacitance [1]. If 
Ej ^ Ec , (j) is well-defined and a phase-coherent Cooper- 
pair current can flow through the junction in the absence 
of an external voltage V. In the opposite limit Ej <C Ec 
an insulating state with a well-defined charge Q on the 
electrodes is possible. At the same time, the dynamics 
of (f) and Q is crucially influenced by dissipation caused 
by the electromagnetic environment surrounding the junc- 
tion. Because of the mutual interplay of quantum mechan- 
ics, nonlinearity and dissipation, the consistent theoretical 
description of Josephson junctions still remains far from 
being complete. 

The dc current-voltage (I-V) characteristics of a 
Josephson junction embedded in a circuit of resistance R 
have been well studied in the so-called overdamped [1] 
case corresponding to small values of R/Rq {Rq = h/Ae^ 
is the resistance quantum) and the ratio Ej/Ec [2-8]. 
For small R < Rq, the supercurrent peak at zero volt- 
age acquires a finite width. With increasing R, quan- 
tum fiuctuations of the phase become more important and 
the supercurrent peak gradually moves to higher voltages. 



This corresponds to the transition (driven by the environ- 
ment) from superconducting behavior found for small R 
to a complete Coulomb blockade when R > Rq. Mean- 
while, the opposite underdamped regime, which has been 
extremely difficult to achieve experimentally, has attracted 
less attention. However, recently experiments were per- 
formed [9,10] on junctions with Ej/Ec > 1, embedded 
in a tunable, highly resistive environment, R ^ Rq, en- 
abling the study of the same junction in different envi- 
ronments. In particular, a voltage peak near zero current 
followed by a back-bending to lower voltages at higher cur- 
rents was observed. This is the so-called Bloch nose [11] 
which, in accordance with a duality property [4,12-14], re- 
sembles the FV characteristic of an overdamped junction 
but with the role of voltage and current interchanged. A 
quantitative comparison between theory and experiment 
has been made in the classical limit where thermal fiuctu- 
ations dominate [10,15]. 

In this Letter we study for the first time the influ- 
ence of quantum fluctuations on the I-V characteristics of 
an underdamped Josephson junction. We show that the 
quantum-mechanical nature of the electromagnetic envi- 
ronment can strongly modify the crossover from Coulomb 
blockade to superconducting behavior. Without fluctua- 
tions the I-V characteristic contains a sharp cusp that con- 
nects two distinct branches: a zero-current finite-voltage 
branch corresponding to Coulomb blockade and a super- 
current branch corresponding to Bloch oscillations of the 
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Fig. 1: I-V characteristics for (7=0.1 and /3)7o=0.1. The dashed 
line corresponds to the classical expression (28) while the dot- 
ted line gives the quasiclassical solution from (29) for /3hu>c='i-- 
Solid lines depict the result from (17) with f3hijJc=W for the 
rounded I-V characteristic and /3hu}c=0-l for the sharp one. 
Inset: (from right to left) the classical limit without fluctua- 
tions (6), the result (28), and the solution from (17) for 5=0.1, 
/3?7o=0.5, and /3!iuic=50. 




Fig. 2: I-V characteristics at T = for huic = 100 Uo- From top 
to bottom, solid lines represent (34) with g=100, 5, 1, 0.5 and 
0.1. Dashed lines are the corresponding curves for the solution 
from (17) with f3Uo=0-5 and f3hLiJc=50. While decreasing R the 
Bloch nose is shifted to the finite current I = $o/2I/, where 
the tunneling of a flux quantum requires an energy exchange 
with the bath. Insets: equivalent circuits of current-biased 
Josephson junction (a) and nanowire-based QPS (b) with the 
diamond symbol representing the phase-slip process. 



voltage. Thermal fluctuations induce transitions between 
these branches thereby smearing the cusp. This ther- 
mal smearing can be described using a Fokkcr-Planck ap- 
proach [15]. Quantum effects show up at temperatures 
comparable with the characteristic cutoff frequency Uc of 
the environmental modes. As we will detail below, at 
intermediate temperatures ksT < fkJc, a quasiclassical 
regime exists. In this regime small quantum corrections 
to the thermal smearing arc found, described by the so- 
called quantum Smoluchowski equation. In the limit of 
low temperatures fc^T <C fkOc, quantum fluctuations not 
only smear the cusp but also shift its position. Specifically, 
with increasing dissipation strength the voltage peak shifts 
to nonzero current: quantum fluctuations induce tunnel- 
ing events of flux quanta that disrupt the Bloch oscilla- 
tions. 

Due to the duality symmetry, our results can be directly 
mapped to the case of an overdamped junction, where 
we have found that the role of quantum fluctuations has 
not been adequately studied, for instance, in the quantum 
Smoluchowski regime [see eq. (29) and discussion below]. 
Finally, we demonstrate that all our results also apply to 
the case of the recently proposed quantum phase-slip junc- 
tions which are realized in superconducting nanowires [16]. 



Model. — Wc start our analysis by considering a 
current-biased Josephson junction (see inset a of fig. 2). 
The junction is shunted by a resistance R and biased by 
an external dc current /;,. As long as quasiparticle excita- 
tions are neglected, the system is described by the follow- 



ing Hamiltonian: 

g = - Ej cos <^ - A J, ^ + H^^ti, , (1) 

where the operators Q and cj) obey the commutation rela- 
tion [(p, Q] = 2ei. The fluctuating charge Qx is associated 
with the current = Qx flowing through the shunting 
resistor. The latter is modeled by a bath of harmonic 
oscillators with frequencies {cj^}, which is described by 
^fbath = {pI + xl) /2. In this model, is rep- 

resented by a weighted sum of oscillator coordinates {xa}, 
XaXa, while the influence of the bath on the 
junction dynamics is entirely determined by the weighted 
spectral function, K{llj) = (7r/2) A^(5(|u;| — uJa). As 
follows from the equations of motion generated by H, the 
choice K{lo) = RQReY{Lu)/2'!rL0 with Y{0) = 1/R repro- 
duces the linear response (Ohm's law) of current to 
the voltage drop V = {h/2e)^ on the resistor, Ix{i^>) = 
Y{LL!)V{iL>) (for Fourier transforms). In the case of interest 
here, Ej » Ec, it is advantageous to switch to the Bloch- 
band description of the Josephson junction [11]. Assum- 
ing that the junction dynamics is conflned to the lowest 
energy band €o{q) = —Uq cos(7rg/e), we arrive at the tight- 
binding (TB) model of our system with the Hamiltonian: 

Htb = -Uo cos{wq/e - x) - {h/2e)Ib(p + i/bath , (2) 

where q is the operator corresponding to the quasicharge 
and Uo = 4^/27^^c (2£j/£c)^^^ exp (^-4:^2Ej/Ec) is 
the Bloch bandwidth. The reduced fluctuating charge 

X = 2n'^XaXa (3) 

a 
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gives rise to the bath correlation function J(r) 
{{x{t) - X(0))x(0)), or explicitly, 



J(t) = 2Rq I 



1 



dw 3T— 



(4) 



where /3 = l/fc^T is the inverse temperature. 

Hamiltonian (2) describes the quasicharge dynamics in 
the lowest Bloch band under the influence of the current 
I = lb — Ix- The group velocity associated with the qua- 
sicharge, deo/dg = (7rt/o/e) sin(7r(j'/e), corresponds to the 
voltage across the junction. Hence for a resistive envi- 
ronment Ohm's law yields Ix = (Vc/ R) sm{Trq/e), where 
Vc = nUo/e is the maximal (critical) voltage the junction 
can sustain. Classically, the problem is equivalent to the 
equation of motion [11], 



q = Ib- {Vc/R) sin(7rg/e) , 



(5) 



describing overdamped quasicharge diffusion with damp- 
ing rate wVc/eR. A stationary solution q = exists if 
h < Vc/R: all the current flows through the resistor, 
while the junction stays in a zero-current Coulomb block- 
ade state with a voltage drop V = Rib- If h > Vc/R, 
q ^ Q and a dynamical state exists at flnite current with 
Bloch oscillations of the voltage. From eq. (5) one finds 
the frequency of the oscillations to be fs = -'^/2e; by di- 
rect integration of eq. (5) over one period one obtains the 
dc voltage [11] 



V = Rh 



(Rhf V^ . 



(6) 



The resulting I-V characteristic corresponds to the afore- 
mentioned Bloch nose and is depicted in fig. 1, inset. 

We now turn to the effect of fluctuations and introduce 
a cutoff of the bath spectrum at frequency Wg, chosen to be 
smaller than the gap ujg = \/2EjEc /h between the lowest 
Bloch bands but higher than the Bloch bandwidth. Fur- 
thermore, the competition between the cutoff energy and 
the temperature determines the nature, rather classical 
or quantum, of the environment. Indeed, starting from a 
classical thermal bath when ksT ^ huc, quantum fluctu- 
ations appear in the quasiclassical region ksT < fkOc and 
become dominant at low temperatures fc^T ^ fkUc- As- 
suming that the effective impedance seen by the junction 
is given by a resistance R in series with an inductance L, 
we can set Y{uj) = 1/{R — iujL), leading to ujc = R/L. 
Then, from evaluating the integral (4) one obtains: 



J(t) = -isign(T) A{t) - M{t) , 



(7) 



with 



A{t) = 7r5(l-e-'^<^l^l) 



(8) 



M(r) = 2g 



J3h 2 \ 2 



|l_e-cM| 



(9) 



where g = Rq/R is the dimensionless conductance, and 
Vk = 2Tik / hp is a Matsubara frequency. 

Before proceeding, we want to comment on the magni- 
tude of the effective inductance L which has been associ- 
ated with the cutoff lJc- The restriction lOc < ujg imposes 
a lower boundary on the values of L to be consistent with 
our assumption of single-band charge dynamics: 



L > 



^/2Ec/Ej 



■'J ) 



(10) 



where Lj = ($o/27r) /Ej is the Josephson inductance 
and $0 = h/2e is the flux quantum. For relatively large 
g > I, the condition (10) can be satisfied even for an au- 
tonomous Josephson junction, where L Lj corresponds 
to the Bloch inductance [17]. In the opposite limit of small 
g <^1, an experimental realization of this model would be 
the use of an environment composed of Josephson junction 
SQUID arrays, whose effective inductance can be tuned to 
relatively large values by magnetic flux [18]. Another op- 
tion includes the use of a quantum phase-slip junction 
where there is no limitation on Wc, as discussed at the end 
of the paper. 

Current-voltage characteristics. — To determine 
the I-V characteristics, we calculate the average of the 
operator for the voltage across the junction 

V{t) = Vc{Uc{t,0)sme{t)), (11) 

where 9 = (ir/e) q — X- The time evolution operator (in 
the interaction picture), 

Ucit, 0) = Tc exp |(iC/o/?i) J^^rJ^s cos0^(t)| , 

(12) 

is deflned along the Keldysh contour with the Keldysh in- 
dex ,s = +/— refering to the forward/backward branch of 
the contour, and Tq denotes the Keldysh time-ordering. 
Expanding Uc in C/q, averaging over the bath and per- 
forming summation over the Keldysh indices yield: 



+ 00 



y(t) = (K/2i)^(-ir(c/o/fi) 



2n+l 



n=0 



/2n+l 



X / dn . . . / dT2„+i V n sinGfc e^" ,(13) 

where for given n the discrete variables fk satisfy 

|/fc+i-M = l (14) 
with the constraint /o = /2n+2 = [5, 19], and 



fe-i 



Gk = Ak,k-lfk + ^2 ^^kk'fk-k' , 

fe'=l 

2n+l 

r„ = -i{n/e)Ib {Tk - Tfc-i) fk 



(15) 



fe=i 
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2n + l 



2n+l fe-1 



^ Mfe,fc_i/2 + 5Z E ^^kk'hh' -(16) 



fc=i 



fe=2 fe' 



Here SMkk' = Mkk' + Mk-i,k'-i - Mfe-i,fe' - Mfe,fe'-i, 
i5^fcfc' = j4fe^fe_fe/_i — ^j. , with the shorthand notation 
Akk' = ^(''■fe — "Tk'), and similar for Mkk'- The above 
equations arc in a form suitable to apply the "nearest- 
neighbor approximation" (NNA), where one assumes that 
SMkk' ~ and SAkk' ~ [5] on the relevant time scale t > 
I/lOc- Then, in the limit t oo, the expression (13) can 
be evaluated exactly, and one obtains for the do voltage V 
across the junction: 

V = Vclm{Wi/Wo) , (17) 

where W„ obey the following recurrent relation (n > 1) 

ia„ (Wn-l - W„+i) = Wn , (18) 



dr sin {nA{T)) e' 



i{v/e)Ii,nT-n M{t) 



(19) 



Eqs. (17)-(19) constitute the central result of this Letter. 
The quantities Wn should be identified with the Fourier 
transform of the quasicharge distribution function W{q) 
in the steady state. 



W{q) = e-'"'^«/^ Wn 



(20) 



with the property W-n = W*. This becomes obvious by 
noticing that eq. (17) can be viewed as a result of aver- 
aging the voltage operator with the quasicharge distribu- 
tion W{q), 



V/Vc = Wq-' J dq W{q) sin (7rg/e) . (21) 

Quasiclassical limit. — To shed light on the range of 
applicability of the NNA for our system, we first study the 
quasiclassical limit where the typical time scale of the qua- 
sicharge dynamics, determined by the damping rate gUo/H 
and the frequency of Bloch oscillations fs [see eq. (5)] , is 
slow. Indeed, if 



gUo/h, fB «: i/np, Wc 



(22) 

one can neglect the terms 0{c^'^''^ , e^"<='^) in eqs. (8) and 
(9). We thus approximate [0(r) is the unit step function] 

A{T)=Trge{\T\-l/uJc), (23) 



and 
with 



M(T) = (ffz.i|r| + A)e(|r|-lK) 

'(3hL0. 



X = 2g 



7- 



27r 



(24) 



(25) 



the approximated A{t) and M(r) is the simplest way to 
provide that they vanish at r = [see the exact expres- 
sions (8) and (9)], which is necessary for the consistency 
of the applied NNA. Assuming also g <^ 1, one obtains 
for a„ entering the recurrent relation (18), 



z e 



-An-' 



an = 



2in — ir] 

where z = (3Uo, V = (e/7r)/3i?4, and [20] 



A = 25 



7 + 



27r 



2tt 



(26) 



(27) 



Inspection of the terms dropped in the NNA reveals that 
in the range (22), where A{t) and M(t) can be simpli- 
fied by their asymptotic expressions, the NNA becomes 
exact in the quasiclassical region where A is small. Note 
that for small g, the NNA can be applicable even at low 
temperature. 

Quantum Smoluchowski equation. — The 

crossover between the classical and quantum limit is 

controlled by the parameter A which is related to quan- 
tum corrections to the position dispersion of a fictitious 
Brownian particle in the harmonic potential [21]. The 
classical limit corresponds to phiOc ^ 1 and g 1, 
or equivalently A, A 0. The recurrences (18) can 
then be solved analytically in terms of the modified 
Besscl functions, Wn = In-irj{z), resulting in the I-V 
characteristics obtained in ref. [15]: 



V = Rh 



sinh (ttt/) 
e/3|/i,(z)|2 



(28) 



Now expanding a„ in eq. (26) to the first order in A 
to include quantum corrections, one can derive from 
eq. (18) a differential equation for the quasicharge dis- 
tribution W{q), 



da 



C{q,d^)-K{e/T,fU"'{q)dl W{q)=Q, (29) 



where U{q) = — [/q cos(7rq/e) — Rhq is a washboard po- 
tential in the charge variable, U'{q) = dqU. The corre- 
sponding Smoluchowski differential operator, 



t{q,d,) = UM+p-^D{q)d,, 



(30) 



is renormalized by quantum fluctuations through both the 
effective potential 

U.n{q) = U{q)+K{e/T:fU"{q) (31) 

and the q-dependent diffusion coefficient 

D{q) = l + 2(3K{elT:fU"{q). (32) 



where 7 = 0.577... is the Euler constant and "(/"(.x) is the Eq. (29) is the so-called Quantum Smoluchowski Equa- 
digamma function. Introducing the short-time cutoff for tion (QSE) [20,21] and constitutes the other main result 
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of tliis Letter. Eq. (29) describes the leading quantum cor- 
rections to the charge dynamics which originate from the 
quantum nature of the bath. According to eq. (27), the in- 
fluence of quantum fluctuations becomes more substantial 
with increasing parameter fSfkOc- The I-V characteristics, 
parametrically dependent on the bias lb and calculated 
for different /3fkOc, are shown in flg. 1. For I3huc ^ 1, 
quantum fluctuations reduce the blockade voltage and fa- 
cilitate a crossover to the Bloch oscillations. Interestingly, 
for small 0huic < 1, the influence of thermal fluctuations 
becomes suppressed by the inertia effect of the "heavy" 
Brownian particle with an effective mass scaled as 
For a wide range of 7b, the junction is locked in the insu- 
lating state because of the lack of energy exchange with 
the environment, which results in a sharp crossover to the 
Bloch oscillations. In other words, reducing the cutoff en- 
ergy below the temperature is equivalent to decoupling the 
junction from the environment (see fig. 1). 

We note in passing that the structure of the QSE ob- 
tained here from the series expansion (17) is different 
from the one derived in ref. [21] for a Brownian parti- 
cle in a slightly anharmonic potential and later applied 
to an overdampcd junction [20]. Therefore we cannot ex- 
ploit the duality property (see below) to treat the Smolu- 
chowski range for underdamped junctions using the results 
of ref. [20] . Another approximation consists of evaluating 
expression (13) using the asymptotes for long times (23) 
and (24) [12]. Such an approach leads to coefficients 
a„ = l3Uo c~-^/^ exp{—2Trgn'^ / (ifkUc) /2i{n — irj) , which tend 
towards the result (26) in the quasiclassical region. As a 
consequence, a QSE equivalent to the solution (29) is also 
recovered. 

Low temperatures. — We proceed by studying 
eqs. (17)-(19) in the low temperature limit, beyond the 
quasiclassical region (22). At very low temperature, 
jShujc 3> 1, the I-V characteristics arc entirely determined 
by the first coefficient ai, eq. (19), consistent with the 
fact that the NNA is exact to the lowest order in Uq. 
Closed-form analytical expressions can be obtained in sev- 
eral cases. For instance, at finite temperatures and for 
small values of g, the I-V characteristic can be written as 

v/v. = u n9 + inpwf ^^^^^^/^^^ ^33^ 

where u = {PUo/4:Tt) {PhcOee^ /2Try^^ and r(x) is the 
Gamma function. The resulting linear resistance R varies 
as a power law with temperature, R cx T'^^s-^) ^ jj^ agree- 
ment with the asymptotic analysis of ref. [4]. At zero- 
temperature one recovers 

V/Vc = (7rC/o/2) P{hh/2e) , (34) 

where P{E) = 1/h J dtexp[J{t) + iEt/h], dual to the 
well-known result for the I-V characteristic for incoher- 
ent Cooper pair tunneling in an overdampcd junction [2]. 
In fig. 2, we plot I-V characteristics for different <?, both at 



zero and at finite temperatures. As g increases, the volt- 
age peak shifts to finite values of the supcrcurrcnt. This 
behavior can be interpreted in terms of incoherent tunnel- 
ing of the phase [22]. Indeed, for small values of g, few 
environmental modes are available. Consequently, only 
elastic tunneling is allowed and a fiat Bloch nose is recov- 
ered. When g is large, the equivalent circuit consists of a 
loop containing the junction closed by the inductance L. 
A phase-slip event occurs when the energy to be released 
^oh corresponds to the energy ^q/2L to add one flux 
quantum $o in the loop, i.e.. when Ii, = ^q/2L. At this 
finite current, phase tunneling disrupts the Bloch oscilla- 
tions and gives rise to a voltage peak. 

Overdcimped Josephson junctions. — A similar 

analysis within the Keldysh formalism can be achieved 
in the case of a voltage-biased overdamped Josephson 
junction in series with a resistance R. This circuit 
is conventionally described by the Hamiltonian H = 
{Q + Qxf /2C - Ej cos (j) - VbQx + -ffbath, where Qx is a 
fluctuating charge on the junction capacitor and Vb is the 
bias voltage. Using an analysis of the equations of mo- 
tion similar to the one preceding eq. (2) and performing 
canonical transformations, one can show that the model 
description of the junction is equivalently given by the 
weak-binding (WB) Hamiltonian 

FwB = -Ej cos (<^ - x/27r) + VbQ + ifbath , (35) 

with X the bath variable deflned as before, eq. (3). In this 

representation, the junction capacitance C is encoded in 
terms of the bath parameters, = {Rq/2tt) ^a^a^ 
while the weighted spectral function of the bath is given by 
KwB{io) ^ 2TrRcZ {lu)/Rquj, with Z-'^{lu) = I/R-iluC. 
The operator for the current flowing through the junction 
is given by / = Ic sin {(p — x/2tt), where Ic = 2eEj/h is 
the critical current. The Hamiltonians Htb and i?wB are 
related by the following transformation: 

(Trg'/e, (/)) ^ {(l>, -nQ/e) , 

(36) 

Uq^Ej, Vb/Rq , X ^ x/27r . 

Consequently, the series expansions for V/Vc in the TB 
model and for I/Ic in the WB model are dual: we can 
transpose our results obtained for the quasicharge dynam- 
ics in an underdamped junction onto the dual case of the 
phase dynamics in an overdampcd junction. A unified ap- 
proach is thus provided, as well as new results such as the 
series expansion for I/Ic (see footnote ^) and the quantum 
Smoluchowski equation. 

Superconducting nanowires. We conclude by 
considering quantum phase-slip (QPS) dynamics in su- 
perconducting nanowires. Based on a duality argument, 

^Thc WB expansion for I/Ic, which is dual to (13), is similar but 
not identical to the one obtained in ref. [7] using a different kind of 
approximation (not the NNA). 
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ref. [16] suggests the following model Hamiltonian to de- 
scribe QPS events: 

Hqps {<t>,Q) = El (0/27r)' - Eg cos (ttQ/c) , (37) 

where Eg is an energy associated with the phase-slip 
process which changes the phase diflFerence cj) over the 
nanowire by 27r, and El = ^l/2L is an inductive energy 
of the wire with a kinetic inductance L. Correspondingly, 
the Hamiltonian of a current-biased QPS junction (see in- 
set b of fig. 2) can be written as 

H = Hqps{<P + Q) + {fi/2e) h 4>^ + i/bath , (38) 

where a fluctuating phase across the junction (j)^ is related 
to the voltage drop Vr = R{Q — lb) on the resistor, (j)^ = 
{2e/h)Vii. It is straightforward to sec that the Hamilto- 
nian (38) is exactly dual to the Hamiltonian of a voltage- 
biased Josephson junction, and can correspondingly be 
mapped onto the Hamiltonian (2) with Uq replaced by Es, 
while the inductance L of nanowire is encoded in terms of 
the bath parameters, = (27r/i?Q) A^Wa- With 
this device, the frequency Wc = R/L results from the 
physical resistance and inductance of the wire, providing 
a natural cutoff of the bath. Our previous analysis for an 
underdamped Josephson junction, and thus the I-V char- 
acteristics, can be directly applied to a superconducting 
nanowire (not being restricted by the TB limit). A typi- 
cal case of nanowire inductance L ~ 1 nH corresponds to 
Wc/27r ~ X 1 THz. Assuming Es/h > 10 GHz [16] 
and T ~ 1 K, as follows from eq. (22), we estimate that 
the QSE range for nanowircs is relevant for g < 0.1. Note 
from the parameters used in fig. 1 that the quantum fluc- 
tuations in QPS jmictions should be substantial. 

Summary. — To summarize, we have studied the in- 
fluence of quantum fluctuations on the I-V characteris- 
tics of an underdamped .Josephson junction. We have ap- 
plied a unifying approach based on the Keldysh formalism 
that enables us to obtain quantitative results for a wide 
range of parameters. Using the NNA approximation, we 
show the significant role of quantum fluctuations revealed 
both in the quasiclassical Smoluchowski regime and in the 
low-temperature quantum regime. In the Smoluchowski 
regime, compared to the case of thermal fluctuations, 
quantum fluctuations mainly lead to a renormalization 
of the parameters describing the quasicharge dynamics of 
the junction [12]. The NNA becomes exact in the case 
of incoherent phase-slip events at low temperatures. In 
this limit, phase tunneling disrupts the Bloch oscillations, 
leading to a voltage peak at finite current. The quantum 
effects are sensitive to both the dissipation strength and 
cutoff frequency, and could be observed in experiments 
as in refs. [9, 10] with a tunable environment. Besides 
Josephson junctions, our results are also relevant for su- 
perconducting nanowircs modeled as quantum phase-slip 
junctions at low temperature. 
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